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Abstract 

Vibrating walls, used to maintain the temperature in a granular gas, modify the system strongly. 
Most conspicuously, the usual one-peak velocity distribution splits into two, asymmetrically po- 
sitioned. A surgical repair of the usual hydrodynamic description is presented that provides an 
account for, and an understanding of, the situation. 
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Loudspeakers operating in air and vibrating walls exciting granular gases are qualitatively 
different devices. The former perturbs the equilibrium state of air slightly, the latter provides 
the largest velocity in the system, and is the very reason why the grains possess kinetic 
energy, why an associated temperature exists. In air at equilibrium, the velocity distribution 
f(v) is a peak of width a/T, centrally located. (We set m, the particle mass, to 1 throughout.) 
A loudspeaker displaces the peak periodically, oscillating it by a velocity much smaller than 
its width, (v) <C y/T. A vibrating wall hits the grains, injecting them into the gas. After 
numerous collisions absorbing much of the initial energy, the grains return to the wall's 
vicinity, to be hit again. Stationarily, when the energy injection equals dissipation, f(v) 
does not change with time. It then has, close to the boundary, two peaks: One for the 
(+) particles leaving the wall swiftly, the other for the (-) particles returning rather more 
slowly. Away from the wall, the two peaks merge into one, because collisions equalize them. 
Yet this circumstance is not a mere boundary problem, as the inefficiency of energy transfer 
limits the size of a granular gas maintained by vibrating walls: If the system is too large, the 
center cools, forms clusters, and is no longer gaseous. If this is to be avoided, the two-peak 
structure extends fairly deep into the gas. That this circumstance is relevant for coming to 
terms with all aspects of granular gases has been emphasized especially by Evesque, see 

The data from • 2D simulations by Herbst et al |2|, who employ boundary conditions 
realistically modeling vibrating walls, and from • micro-gravity experiments obtained in an 
Novespace Airbus (Campaign 2006) in a 2D vibro-fluidized granular system ^j], show two 
intriguing results that cry out to be understood. First, instead of being characterized by a 
pressure as any gas, the system, seemingly emulating a solid, develops an anisotropic stress - 
such that o xx normal to the pair of vibrating walls is strictly constant, while the tangential 



component o yy [x) depends on x - though force balance is satisfied [2]. (a xx is constant 
except in a genuine boundary layer of 1.5 grain diameters.) Second, the temperatures T x 
and T y , for the velocity distribution along x and y, are different. The reason is probably 
that it takes many lossy collisions before the injected momentum along x gets fed into the 
kinetic energy along y, or z. But the questions remain where the solid-like behavior comes 
from, and how to establish an equation of state that separately relates a xx and <J yy (x) to the 
temperature and density fields, reproducing especially a xx s highly accurate constancy. 

Conventional hydrodynamics employs as variables the densities of mass p and momentum 
pv, both conserved, in addition to T, a measure for the averaged kinetic energy per particle 



in dilute systems, and more generally, for the averaged total energy in microscopic degrees 
of freedom. The pressure is given as P = pT in an ideal gas, and more generally as p = 



(pl-l)( w -Ts) 



where w, s are respectively the energy and entropy density. Granular 



hydrodynamics 3] , as first proposed by Haff [7], has the same variables, though T relaxes, 
going to zero if the grains' kinetic energy is not replenished. 

To account for a two-peak fluid, it seems obvious that one should employ as variables two 
different sets of T^p, pv, for the (+) and (-) particles, along with formulas for the pressures, 
P + and P_, see [8y. In addition, one would need two additional sets of T,p,pv for y and 
z. We are afraid this opens a Pandora box of state variables, neither conserved nor truly 
independent, rendering the resultant theory arbitrary and unwieldy. One should instead, we 
believe, retain the conserved variables p and pv, also T - as the average width of all peaks 
in the system jjj. It should then suffice, for a minimal, surgical modification, to introduce 
two additional variables. 

The first is the distance between the peaks, which is the crucial second length of a two- 
peak distribution, much more relevant than the difference between the two widths. (In a 
single-peak distribution, the width is the only length scale. But there are many more in 
a two-peak one, and providing the two widths alone is not sufficient to characterize the 
distribution.) Being a velocity difference, the new variable Aj is odd under time reversal 
and a vector. In our case, only A x ^ 0. 

Second is the difference between T x and T y , between the average width along x and the 
only width along y. More generally, we have 5Ti = Ti — T, i — x,y, z, with £V 5Ti = 0. 
These are - similar to the order parameter of nematic liquid crystals [10(] - the diagonal 
elements of a symmetric, traceless tensor. (They do not form a vector, because T x does not 
distinguish between x and —x.) So a tensor, even under time reversal, needs to be added. 
Instead of STij, however, we employ t^, the deviation of the granular temperature T g as 



considered below, see Eqs ( fTTIoT) . In our case, on 



Granular solid hydrodynamics (GSH) Hi, LL3J was derived employing the hydrodynamic 



y txx tyy 7^ ^' 



procedure. Relying on general principles valid irrespective how dense or rarefied the system 
is, it leads to equations that include collisions and enduring contacts, are valid both in the 
dense, elasto-plastic limit, and the rarefied one. For instance, GSH is capable of accounting 
for the relaxation of the temperature until it is zero. 

Although the present system deviates from rarefied gas in the opposite direction, towards 



ballisticity, the hydro dynamic procedure still works, if we add variables that characterize 
the deviation from local equilibrium. The reason is the hierarchy of equilibria: Although 
the two peaks or three widths are not in equilibrium with one another, the elements within 
each are well thermalized. 

In deriving GSH, a granular heat w is introduced - with s g the granular entropy, and T g = 
dw/dsg the associated temperature. It quantifies the energy contained in the mesoscopic, 
intergranular degrees of freedom, especially the strongly fluctuating part of the grains' kinetic 
and elastic energy. Expanding w in s g , we have 

w = s 2 J(2bp) = bpT 2 /2, b~(p cp -p) a \ (1) 

with ai = const. The lowest order term is quadratic because equilibrium, or minimal energy 
w = 0, is given for s g = bpT g = 0. (This is quite the same idea as with any Ginzburg-Landau 
energy functional, just without the fourth order term, or a phase transition.) The density 
dependence of b(p), with p cp the random close density, is chosen such that the associated 



pressure 
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li, P = {p-l -l)(w- T g s g ) = -\p 2 T 2 db/dp, given as 



'dp ^v" 7 "rs/ 2^ 's' 

_ axpw _ \a\p 2 bT 2 

Pep - P Pep - P 

is appropriate for all densities, see eg. 



(2) 



14] 



For a rarefied system, in which the elastic contribution to the energy is negligible, we 
may identifying the energy w = \bpT 2 with the kinetic energy per unit volume pT (for a 2D 
system), implying bT 2 ~ 2T, especially in the above expression for the pressure. Clearly, 
taking T g ~ \/T, Haff's granular hydrodynamics is retrieved. 

If the temperature is maintained by vibrating walls, we need (as discussed above) A, and 
tij as additional variables. They also contribute to the energy which, in an expansion in all 
three variables, becomes 

w = (bpT 2 + c P A 2 + eptl)/2. (3) 

These variables relax, specifically because they possess energy that may be redistributed 
among microscopic, inner-granular degrees of freedom (such as phonons). The energy being 
quadratic, the relaxation stops when the variables are zero, and the energy vanishes. 

Taking c, e = const, independent of the density, the pressure is not changed by introducing 
the new variables, and remains as given in Eq (T5]). (Same as with P ~ db/dp, additional 



pressure contributions would have resulted from dc/dp and de/dp.) We assume this for 
simplicity, as we are more interested in an anisotropic stress, less in modifying a given 
pressure. 

Next we relate A x , t xx to parameters of f(v x ) and f(v y ), the velocity distributions, as the 
latter are independently measurable. Denoting the norm as N = y/TT%(l + a) and 2k B = 1, 
we take 



f(v x ) 



a exp — h exp 



(4) 



and f(v y ) = f(v x -+v y ,T x ->Ty,£ = 0), with (v x ), (v y ) = 0, see Fig[U The energies along x,y 
are then w x = \pT x + pa£ 2 and w,, = \pT y) implying, first of all, a— >c, £— ^A x . 

Taking = |(T + ST XX ), T y = h(T + ST yy ), with 5T XX + 5T yy = 0, £ = 0, we find that 
the energy, Wj, + ^ = |p7", does not depend on But taking 



T x = b(T g + t xx f/A, T y = b(T g - t xx f/A, 



(5) 



we obtain w x + w y = ^pb(Tg + t xx ), as in Eq (j3J). This discrepancy may be surprising at 
first, but results from pT being the energy of an ideal gas, or the kinetic energy of a rarefied 
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FIG. 1. Velocity distribution as simulated in Q|, measured in and parameterized by Eq @. 
Fig (a) and (b) show f(v x ), for the velocity perpendicular to the vibrating walls, with (a) showing 
/ close to the wall, and (b) showing / in the middle. Fig (c) shows f(v y ), while Fig (d) again 
shows f(v x ). Symbols are measurements, dotted lines simulation, and full curves are Eq (T4|). For 
the two-peak distributions, we have T x = 10 and 1, £ = —2.2 and -1,37, a = 2.3 and 3, for Fig (a) 
and (d), respectively. 
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gas, with no collisional contributions. Yet collisions are what equalize T x and T y . On the 
other hand, w = \pbT^ is the total energy of an interacting system. Increasing T g by t xx in 
one population, and decreasing it by t xx in another, must lead to an energy increase, as t xx 
would not relax otherwise. Still, we should not take the equality of b and e seriously, as it 
hinges on the precarious assumption that w x , w y remain sensible quantities in an interacting 
system. In contrast, expanding w in tij leading to Eq ([3]) is a generally valid approach. 
Having specified the additional variables and their contributions to the energy, we may 



employ the same hydrodynamic procedure as used for GSH |l3j, to set up their equations 
of motion. The generalized hydrodynamics consists of continuity equations for momentum, 
dt(pvi) + V^crjj + pViVj) = 0, and mass, d t p + Vj(/w,) = 0, in addition to the balance 
equations, 

T g\ d tSg + Vi(s g Vi - KgS/iTg)} = VgVtjV-j - l g T^ (6) 
d t Ai + Vj(AiVj - K A VjAi) = aAjvfj - 7aA;, (7) 
dtUj + V k (UjV k - KtVkUj) = - jtUj, (8) 
with cry = P5ij — rjv*j — cp a AjA^ — epfttij. (9) 

The first equation is the same as in GSH. It sports a convective, a diffusive (~ n g ) and a 
relaxative term (~ j g ), in addition to viscous heating, with r\ g the viscosity, and v*j the 
shear rate - v*j being the traceless part of = |(Vjfj + Vjfj). Employing T g ~ yT, one 



sees that Eq (Q is the same as Haff's energy balance [151 ] . Eqs fl71l8|) are new, but quite 
similar to (Q. They also each sport a convective, diffusive, and relaxative term. Instead 
of viscous heating, however, there is a linear, off-diagonal Onsager term: with v*j as the 
thermodynamic force, Aj as the preferred direction, and a an Onsager coefficient in Eq (J7j); 
with v*j as the force, no preferred direction, and /3 another Onsager coefficient in Eq (jSJ). 

The stress cr^ consists of pressure, viscous stress (with bulk viscosity neglected), and 
the two counter Onsager terms. The signs of a, (3 in the three equations obey Onsager 
reciprocity relation; and because both AjA^u^- and are odd under time inversion, their 
respective contribution to the production of true entropy (not displayed) vanish. There is 
no constraint on the sign or magnitude of a, (3; both are functions of the density. 

To solve Eqs ( 161171181191 . we first note that for the above discussed Herbst geometry, setting 
Vi,Vij = 0, and assuming dependence only along x, we have a xy = 0, 

a xx = P - cpaA 2 x - ep(3t xx , a yy = P(x) + ep(3t xx . (10) 



Force balance VjO"jj = requires a xx = const, but leaves a yy undetermined. Denoting 



n2 — «9 p2 — n A p2 



A = l\ = with 2L the distance between the two vibrating walls, and 



9 — 2j g ' A 2 7A ' t — 2-y. 



employing the boundary conditions: T g = T , A x = A , t xx = t at x = 0, and T 9 = T , 
A x . = — A , t xx = t at x = 2L, the stationary solution, for d t s g , d t Ai, d t Uj = 0, is 



T g _ coshj 7T *** _ cosh^ A x _ sinh^ 
T cosh ^ ' to cosh ^- ' A sinh ^ 



(11) 



Note we have taken all transport coefficients, generally functions of p, T g , Aj, as constant. 
(Although k s ,7 s ~ T g , see [Hi, [l5j, this does not change the solution if included, since the 
equation contains only the ratio l 2 g = ^-.) Searching for an understanding at present, we 
are unabashedly qualitative. Given the scarcity of experimental and simulation data, there 
is too much arbitrariness for more quantitative considerations. For a comparison of theory, 
experiment, and simulation, see Fig [2|3l 

Summary: The physics of a granular gas has some idiosyncratic features, if its tempera- 
ture is maintained by vibrating walls. These are most notably first a two-peak distribution 
for the velocity perpendicular to the vibrating walls, second a much narrower peak width for 
the velocity along the walls, and third, as a result of the above two, a solid-like anisotropic 
stress. We modified the hydrodynamic theory by introducing two additional variables: A x 
for the distance between the two peaks, and t xx = —t yy for the difference between the widths, 
where the former is part of a vector, and the latter part of a symmetric, traceless tensor. Sim- 
ilarly to the granular temperature T g , both variables characterize the velocity distribution 
of an inelastic gas. All three diffuse and relax, displaying a similar macroscopic behavior. 
While the scalar T g gives rise to an hydrodynamic pressure, the vector Aj and the tensor 
tij contribute to an anisotropic stress. Even under strongly simplifying assumptions, the 
calculated stress displays remarkable resemblance to that of simulations [2] and micrograv- 
ity experiments 3|. We conclude that for this case, the basic approach of a hydrodynamic 
description remains valid. 
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FIG. 2. Variation of stress and kinetic energy along x. We employ a xx = const as the unit of stress 
and energy density. Full curves are hydrodynamic results, taking w x = pbTg / '4 + pet xx /2 + pcA x /2, 



and w y = pbT'/A + pet^/2, with b ,c,e = 1, and i g = 1,1.1,2,1, t t = 0.8,0.9,1,0.5, £ A = 
0.18,0.28,0.55,0.3, m = 1.38,1.2,0.33,0.9, p cp = 0.7,0.7,0.7,0.8745, -a = 13,1.49,0.772,0.5, 
— f3 = 3.68,2.14,2.18,3.5 for (a), (b), (c), (d), respectively. Symbols are from micro-gravity mea- 
surements of [3J, and dotted lines from simulations of 2|. 

9 




10 



